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ABSTRACT 

This work presents high S/N spectroscopic observations of a representative sample of nearby ellipti- 
cal galaxies. These observations provide a strong test of models for the formation of elliptical galaxies 
and their star formation histories. Combining these data with the Gonzalez (1993) data set, a volume 
limited sample of 45 galaxies has been defined. Results are in agreement with previous studies: the 
existence of the metallicity hyper-plane and the Z-plane of Trager et al. (2000) is confirmed, and the 
distribution is clearly due to physical variations in stellar population parameters and not measurement 
uncertainty. Trends between stellar population parameters and galaxy structural parameters suggest 
that angular momentum may determine the chemical abundance of a galaxy at a given mass. 
Subject headings: galaxies: elliptical and lenticular, cD — galaxies: abundances — galaxies: stellar 
content — galaxies: formation — galaxies: evolution — galaxies: general 



1. INTRODUCTION 

Onc e thought to be simp l e old and metal-rich sys- 
tems (Ba ade fe GaposchkirJ I1963D . elliptical galaxies 
have been studied extensively over the past thirty years. 
A number of studies (e. g. Gonzalez 1993, hereafter G93; 
Rose et al. 1994; Jorgensen 1997; Tantalo, Chiosi, & 
Bressan 1998; Kuntschncr & Davies 1998; Trager et al. 
2000a; Longhetti et al. 2000; Terlevich & Forbes 2002; 
Caldwell, Rose, & Concannon 2003) have shown that el- 
liptical galaxies span a wide range of ages, with a rel- 
atively small spread in metallicity [Z/H]. The enhance- 
ment ratio, [a/Fe], of Mg and nucleosynthetically related 
elements relative to iron peak elements is also found to 
b e significantly greater than solar. 

ITrager et all l|2000bl hereafter TFWG) found that el- 
liptical galaxies occupy a two-dimensional "metallicity 
hyper-plane" in the space of [Z/H], [a/Fe], log t, and 
log tr, where t is the light-weighted single stellar pop- 
ulation (SSP) equivalent age and a is the velocity dis- 
persion. Although cluster galaxies are older on average 
than field gala xies, a wide rang e of ages is seen in both 
environments l).IorgenserJll999[f . iGebhardt et alJ l|2003f) 
suggested that a frosting of continuing star formation is 
necessary to explain the lack of color evolution in early 
type field galaxies since z = 1. TFWG also proposed 
a frosting model of star formation, as will be discussed 
in more detail in § 3.4. This result emphasizes the im- 
portant point that the analysis of galaxy properties us- 
ing SSP models does not accurately reflect the actual 
star formation historie s of most galaxies. As TFWG and 
IGebhardt et al.l l)2003j) also note, the observables do not 
discriminate well between different star formation his- 
tories. In multiple stellar population models, different 
combinations of mass fractions, ages, and metallicities 
can produce observationally indistinguishable galaxies. 

Because the a elements are primarily produced in SN 
II while SN la produce mostly Fe, [a/Fe] may reflect the 
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duration of star formation assuming a constant, universal 
IMF. In high-redshift starburst, [a/Fe] will be high, as no 
SN la could have gone off to introduce additional Fe. In 
contrast, the long and continuous star formation histories 
typical of spiral galaxies produce [a /Fe] ~ as a result 
of plentiful SN la enrichment over billions of years. Thus 
abundance ratios might constrain the formation mech- 
anism of individual galaxies. Two extreme formation 
mechanisms are a fast clumpy collapse model and a spi- 
ral merger model. The fast clump y collapse mod el is in 
essence a refinement of the classical lLarso mono- 
lithic collapse scenario: protogalactic gas clouds collide 
and merge at high redshift, forming stars rapidly as the 
gas dissipatively collapses to the center of the system. In 
a major merger of two spirals, the bulk of the total mass 
is already in stars formed over long timescales, while the 
remaining gas is metal- rich with near solar [a/Fe]. An 
elliptical formed in a fast clumpy collapse would have 
high global [a/Fe] values since most stars formed from 
gas enriched by only a handful of SN II. 

Previous spectroscopic data sets (e. g., G93, 
iKuntschner fe Daviesl lT998') do not necessarily include a 
representative sample of elliptical galaxies in the local 
universe. Thus conclusions based on such data sets may 
reflect some unanticipated bias in sample selection rather 
than physical relations between structural and/or stel- 
lar population parameters of elliptical galaxies. Stellar 
population modeling has advanced considerably in recent 
years: the inclusion of nonso lar abundance ratios a mong 
q-elements a s pionee red by iTripicco fc Belli (|1995l) and 
ITrager et alJ ll2000al) has been included in more re- 
cent models such as lThomas. Maraston. fe Benderi ( 2003, 
hereafter TMB). 

This paper is organized as follows. Sample selection, 
observations, and data reduction are described in § 2. 
The data are analyzed and stellar population parameters 
are derived in § 3. Conclusions are presented in § 4. 

2. THE DATA 
2.1. Sample Selection 
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A volume-limited sample is defined, consisting of all 
45 early type (E or SO) galaxies with m — M < 32.5 (v < 
2200 km a- 1 for H = 70 km s" 1 Mpc -1 ), M B < -19.5, 
5 > —20°, and \b\ > 15°. The galaxies were selected 
using the NASA/IPAC Extragalactic Database (NED). 
The cumulative luminosity distribution of this sample is 
consi stent with that of a Schechter function (Schechter 
1976) with the standard parameters of M* = —20.5 and 
a = —1.07. This suggests that the sample proposed here 
provides a reasonable sampling of nearby elliptical galax- 
ies within a magnitude of M* , and that the sample does 
not suffer from significant incompleteness near the faint 
limit. The dist ance measurements used in selecting this 
sample are from lTonrv et all 1)200 If) . The volume-limited 
sample is presented in Tabled Many galaxies were pre- 
viously observed by G93, and are noted as such. For 
newly observed galaxies we list S/N per pixel near the 
H/3 line. Although the sample has somewhat lower S/N 
than the G93 sample, this has only a small effect on the 
total observational errors due to the typical size of cali- 
bration uncertainties; see §2.6. The 10 galaxies that were 
not observed are consistent with being drawn from the 
same population as the 35 galaxies that were observed, 
as judged by one-dimensional Kolmogorov-Smirnov (K- 
S) tests on their distribution in velocity dispersion and 
Mb magnitude. 

2.2. Observations 

The observing program used the KAST spectrograph 
with the 1200 lines/mm grating blazed at 5000 A on the 
Lick 3 meter telescope. Observations consist of 4 — 6 
25 minute exposures on a galaxy, usually with the 145"- 
long slit oriented along the galaxy's major axis, inter- 
spersed with 5 minute "blank" sky exposures on fields 
several arcminutes away from the target galaxy. The 
spectral range 4200A-5600A was observed, with instru- 
mental resolution of approximately 100 km/s. The veloc- 
ity dispersions of the galaxies are typically ~ 200 km/s. 
Early observing runs suffered from slightly poorer reso- 
lution as a result of suboptimal spectrograph focus. The 
slit width was typically 1.5", though in poor seeing con- 
ditions this was increased to 2". Spectral resolution was 
measured independently for each observing run to ac- 
count for changes in resolution due to different slit widths 
or focus settings. The plate scale for the spectrograph 
in this configuration was 1.17A/pixel in the dispersion 
direction and 0.8"/pixel in the spatial direction. 

Numerous Lick index standard star observations were 
made, drawing from the s et of stars in common between 
the Lick/IDS standards of lWorthev et all <|1994H and the 
I.TonesI l)1999f) standards. 

2.3. Data Reduction 

Flat fields were created for each night of each observ- 
ing run. The dome flat input images were averaged in 
the spatial direction to one dimension, and this median 
dispersion axis was then expanded back to the full size 
of the CCD. The final flat field images used in this work 
are the input flat fields divided by the one dimensional 
median for each row of the CCD. These flat field images 
correct for pixel-to-pixel quantum efficiency variations 
only. The data do not extend to red wavelengths where 
CCD fringing corrections become important. 



Proper sky subtraction is crucial for accurately mea- 
suring line strengths in galaxies. Oversubtraction of sky 
flux creates spuriously large line strengths, while under- 
subtraction has the reverse effect. For all but the most 
extended galaxies, the ends of the slit provide adequate 
sky subtraction. This simple technique was used where 
applicable for the galaxies observed during or after 2002. 
The most extended galaxies, and those observed prior 
to 2002, were reduced using the offset sky spectra in- 
terspersed with the galaxy observations. The observa- 
tions were not taken under photometric conditions, and 
the sky was found to vary unpredictably over timescales 
much shorter than that of a typical galaxy exposure. 
As a result, averaging the sky spectra taken immedi- 
ately before and after each galaxy spectrum and scal- 
ing the sky spectrum by the ratio of exposure times is 
not sufficient for accurate sky subtraction. The proper 
amount of sky continuum was estimated by forcing the 
sky-subtracted radial b rightness profile to c onform as 
closely to an r 1 / 4 law m^^yc^lgma ^948) as possi- 
ble, as explained below. iBurkertl l)1993f) showed that the 
r 1 / 4 law provides an excellent fit to all elliptical galaxies 
between 0.1 < r/r e < 1.5. 

First, the IRAF 1 task GEOMAP was used to rectify the 
galaxy spectrum such that a fixed spatial location along 
the slit lies at a constant y pixel coordinate. Cosmic rays 
and bad pixels were removed at this step. The median of 
the galaxy spectrum images was subtracted from each in- 
dividual frame. The resulting difference frames were then 
divided by the median image plus a floor value equal to 
ten times the standard deviation at the edges of the slit. 
The resulting fractional deviation image was used to flag 
deviant pixels by interactively setting a threshold value. 
The flagged pixels and all neighboring pixels were masked 
from the input images, which were then combined to 
form the galaxy spectrum. Each input exposure was cor- 
rected for frame to frame variations in the background 
flux, putting every exposure on a common relative sky 
continuum scale. Specifically, the combined galaxy spec- 
trum was subtracted from each cosmic ray corrected in- 
put spectrum. The statistics of the difference image near 
each edge of the slit were used to estimate a relative 
sky offset between the input images. The measured off- 
sets were added or subtracted from each spectrum to 
place them all on a common relative sky continuum scale 
prior to measuring the absolute sky subtraction amount 
needed to fit an r 1 / 4 profile. Since sky intensity varies 
slowly with wavelength, a wavelength-independent offset 
was used in order to best match the sky continuum lev- 
els near the H/3, Mg&, Fe5270, and Fe5335 lines. Next, 
the combined galaxy spectrum was collapsed along the 
dispersion axis in two wavelength ranges, one around the 
wavelength of H/3 and one around the wavelengths of the 
Mgb, Fe5270, and Fe5335 lines, to create radial intensity 
profiles. An r 1 / 4 profile was fit to the galaxy radial pro- 
files at intermediate radii. Finally, a constant offset was 
subtracted across the entire spectrum. This sky contin- 
uum level was determined to be that offset which results 
in the best fit to a de Vaucouleurs profile extending to 

1 IRAF is distributed by the National Optical Astronomy Ob- 
servatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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large radii. In practice little difference was found be- 
tween the best fit offset in the two wavelength ranges 
used to produce radial intensity profiles. 

Extraction and dispersion correction of the spectra 
used standard IRAF packages. The input spectra to 
these tasks incorporate the relative and absolute sky con- 
tinuum adjustments described above, and have been me- 
dian combined with cosmic ray rejection as above, but 
have not been rectified. Instead, the extraction proce- 
dure traces the object spectrum across the CCD. The 
galaxy spectra were extracted in apertures extending to 
r e /8 from the g alaxy center as per G93. Effective radii 
were taken from lFaber et all l)1989fl . Galaxies were dis- 
persion corrected using lamp spectra taken immediately 
before or after the galaxy observations, at the same tele- 
scope position. These lamp spectra were extracted using 
the same trace as the galaxy spectra. A high-order poly- 
nomial (typically fifth order) was fit to the ~ 30 lamp 
lines, and this dispersion correction was then applied to 
the galaxy spectrum. Typical wavelength calibration un- 
certainties are ~ O.IOA. 

2.4. Lick Indices 

Line wid ths of the Lick indices defined in lTrager et alJ 
(|1998J) and lWorthev fc Ottavianl l|1997|) were calculated 
using a ver sion of th e bwid program, provided by 
R. M. Rich l|Ricblll998fl . As bwid is not an IRAF pro- 
cedure, the data were read out into ascii files for input. 
Test spectra from Guy Worthey's web page 2 were used to 
confirm the accuracy of bwid. Two sets of Lick/IDS in- 
dices were measured. The primary set of indices consists 
of H/3, Mg6, Fe5270, and Fe5335, while H 7f , Ca4227, 
G4300, Fe4383, Ca4455, Fe4531, C 2 4668, Fe5015, and 
Fe5406 form the secondary set of indices (Table EJ- The 
two sets are distinguished by both scientific and practical 
considerations. The primary set comprises the indices 
in common use in previous observational studies (e. g. 
iTrager et al.ll2000atl. Fur t hermo re, the spectral range of 
the models of Sch iavonl l|2005|) allow for sophisticated 
error simulations of the primary indices, while simpler 
error estimates must be used for many of the secondary 
indices; see § 2.5 for details. 

Several corrections must be applied before equivalent 
widths can be measured. First, the galaxy sample must 
be adjusted to a comm on radial velocity , v = 0. The ex- 
tensive model library of lSchiavonl l|2005|) was used to pro- 
vide an empirical Doppler correction. The galaxy spectra 
were cross correlated with an old, metal rich model spec- 
trum, and the resulting velocity was used in the Doppler 
correction. 

Next, the galaxy sample was corrected to zero ve- 
locity dispersion. The following procedure was used. 
The model spectra were rebinned to match the Doppler- 
corrected galaxy spectrum. An analytic continuum func- 
tion was fit to the galaxy and each model spectrum. The 
model spectra were divided by their respective contin- 
uum functions, then multiplied by the galaxy continuum 
function to match the shape of the observed galaxy spec- 
trum. A copy of each model spectrum was then smoothed 
to the velocity dispersion of t he gal axy using the a mea- 
surements from iFaber et alJ l)1989|) . A least-squares fit 
was performed to select the appropriate model for each 

2 http://astro.wsu.edu/worthey/html/system.html 



galaxy. The fit was performed within a spectral region 
surrounding the Mg b, Fe5270, and Fe5335 lines: 5140- 
5370A. To empirically refine the galaxy's velocity dis- 
persion, the a = version of the best fitting model 
spect rum was used as a template for the Pixfit code 
from Ivan der Marell l)1994[) . The resulting a was then 
corrected for the different instrumental resolution of the 
galaxy and model spectra. The best fit model was then 
rebinned to a logarithmic wavelength scale and gaus- 
sian smoothed to the measured velocity dispersion of the 
galaxy, again taking the intrinsic resolution of the model 
spectrum into account. After transforming back to a lin- 
ear wavelength scale, line widths were measured in the 
a = model spectrum smoothed to Lick/IDS resolution 
and the model spectrum smoothed to the galaxy's veloc- 
ity dispersion and then to the Lick/IDS resolution. The 
ratio of the two measurements is the correction factor 
applied to the corresponding equivalent width measured 
from the galaxy spectrum. For a few indices (Fe4383, 
Ca4455, Fe4531, and C24668) no model spectra were 
available covering the appropriate waveleng ths; for these 
indice s cr-corrections were calculated from ITrager et alJ 
l(T998h . 

As indicated above, it is necessary to degrade spec- 
tra to the resolution of the Lick Image Dissector Scan- 
ner instrument which was used to make the stellar ob- 
servations upon which the Lick/IDS system is based. 
The resolution function of th is instrument is taken from 
iWorthev fc Ottavianil l)1997|) . Since the resolution varies 
strongly with wavelength at the blue end, smoothing was 
done piecewise in 100A increments in that part of the 
spectrum, while a constant value was appropriate for the 
red end. 

An emission correction was required to accurately mea- 
sure the H/3, H 7F , and Fe5015 indices. The [OIII]A5007 
line was used to estimate the amount of H/3 emission 
(Trager et al. 2000a). Since [OIII]A5007 is adjacent to 
the FeA5015 line, it is impossible to accurately deter- 
mine the continuum level. Instead, model spectra were 
used. The best fitting model was subtracted from the 
galaxy spectrum, and the intensity of [OIII]A5007 was 
measured on the difference spectrum. This line was mea- 
sured by hand using the IRAF task SPLOT. The re- 
sult was then converted to equivalent width using an 
estimate of the intensity level of the galaxy spectrum 
near [OIIIJA5007. The [OIII] emission w as converted 
to an H /3 correction using Equation 2 of ITrager et alJ 
(2000g j). H^F and Fe5015 were corrected as described 
in iKuntschner et alJ ((2002), with correction factors of 
0.36 x [OIIL1A5007 and 0.61 x [OIII]A5007 respectively. 
In NGC 1052 H/3 appears only in emission. NGC 3226, 
with the next strongest [OIIIJA5007 line, also shows H/3 
emission within the aborption line, while all other galax- 
ies show at most a slight distortion in the shape of the 
H/3 absorption. 

Lick index standard stars were used to calibrate the 
output line widths. The reduction procedure for stan- 
dard star observations is a simplified version of that de- 
scribed above for the galaxy observations. A typical stan- 
dard star is observed in five consecutive exposures of a 
few seconds integration time. After flat fielding and dark 
subtraction, the individual frames are median combined 
and then extracted. Sky subtraction is performed via a 
linear fit between two sky regions near the star on the 
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slit. The stellar spectrum is extracted in a user-defined 
aperture extending to the points at which the PSF dis- 
appears into the sky background. This aperture is traced 
across the CCD as above. As with the galaxy spectra, 
the accompanying calibration lamp spectra are extracted 
in the exact same aperture as the stellar spectrum. After 
dispersion calibration, the smoothing described above is 
applied. Radial velocities are determined empirically by 
cross correlation, in this case using the Jones spectra of 
the same star as templates. The standard star spectrum 
is read out into an ascii file, and line widths are measured 
by bwid. Corrections for each of the Lick indices are de- 
rived by comparing the observed stellar line widths for 
each observing run to the standard values. These calibra- 
tion offsets are not constant from run to run, particularly 
for the Fe5335 index. 

The data from this project are not flux-calibrated, 
while the Lick/IDS system was calibrated to a quartz 
lamp response curve. This mismatch is absorbed into the 
calibration corrections. The standard stars and galaxy 
spectra were observed with the same CCD, grating and 
grating tilt, and are thus on a common system. Errors in 
the mean for the IDS calibration to the primary indices 
are 0.04-0.06A, though run to run variations are often 
much larger. Fully corrected measurements for each in- 
dex, including estimates of the total error (§ 2.7) in each 
index, are presented in Table [3 

Figure n shows the measurements for each index plot- 
ted against the galaxy's veloc ity dispersion. The Mg- 
o~ relation i|Bender et alJll998Tl is clearly evident, with 
Ca4227 and C24668 following the same trend. The iron 
indices show similar behavior but with greater scatter, 
forming a mass-metallicity relation. The two Balmer in- 
dices both follow the well-known trend that larger galax- 
ies tend to have smaller values than smaller galaxies. The 
five galaxies with large (> 0.5A) [OIII]A5007 emission 
are omitted from subsequent analysis due to their large 
uncertainties. 

2.5. Error Simulations 

Empirical error simulations were performed as follows 
using the model spectra described previously. The model 
spectrum used in the velocity dispersion correction forms 
the basis for the error simulations for the corresponding 
galaxy. 

Several potential sources of random error were simu- 
lated. Poisson noise was modeled by adding the quadra- 
ture sum of the galaxy poisson noise and sky spectrum 
poisson noise to the model spectrum. The poisson errors 
were derived from the raw spectra (flat fielded only, no 
dark or sky subtraction), in units of electrons rather than 
ADU. 

Error in the subtraction of the night sky lines was es- 
timated as follows. The pixel value of the sky lines was 
multiplied by a small factor / and added in quadrature 
to the above poisson error. The value of / was chosen 
such that the simulated spectrum best matched the true 
galaxy spectrum. Sky line subtraction error proved to be 
unimportant: / = provides a good match to the galaxy 
spectra. 

To investigate the possibility of flat fielding error, a 
noise floor was introduced into the error calculations. 
This floor was simply added in quadrature to the poisson 
error, and the value of the noise floor was chosen so that 



the simulated spectrum best matched the galaxy spec- 
trum. As with sky line subtraction noise, flat fielding 
noise proved to be unimportant, and a noise floor of zero 
was used. 

Two other sources of error are directly measured and 
included in the simulations. The uncertainty in the dis- 
persion solution is calculated by IRAF, and the error in 
sky continuum subtraction is easily estimated. These 
errors are included in the simulations, the first as a ran- 
dom shift in the wavelength scale and the second as a 
random shift in the continuum level of the model spec- 
trum. Sky subtraction error may be underestimated by 
this method, particularly at the blue end of the spec- 
trum, far from the wavelength range in which the sky 
subtraction measurements were made. 

Error simulations were performed in groups of 1000 to 
provide accurate statistics. Each simulation spectrum 
was measured by bwid in exactly the same way that ac- 
tual galaxy spectra were measured. Three galaxies were 
used as the basis for these simulations: NGC 3115, the 
galaxy with the highest S/N in the sample, NGC 1209, a 
galaxy with S/N typical for the bulk of the galaxy sam- 
ple, and NGC 4168, a galaxy with relatively poor S/N. 

2.6. Calibration Uncertainty 

The uncertainty in the calibration to the Lick/IDS 
system deserves particular attention because for many 
indices this is the dominant source of error. The ob- 
served standard star primary index measurements are 
compared to the Lick/IDS values in Figs. [3 and The 
former shows that a constant calibration offset can cor- 
rect the observations onto the standard system — no 
first-order term in index strength is needed for the cal- 
ibration. Fig. ^ highlights the run-to-run variations in 
the necessary calibration offset. These offsets are calcu- 
lated separately for each run, except for 2003 January 
2 (Run 4) which was incorporated in the 2002 Novem- 
ber (Run 3) run calibration due to the limited data ob- 
tained that night. The calibrated standard star index 
measurements are compared to the Lick/IDS values in 
Figs. El and plotting the difference between calibrated 
and standard values against the standard value and the 
run number, respectively. The external error in calibrat- 
ing the data to the Lick/IDS system is taken to be the 
error in the mean calibration offset for each run. These 
external errors are shown at the right of Figs. 01 and |3J 
the error bars are omitted in Figs. and for clarity. 

The repeated observations of a subset of standard stars 
in several observing runs allows the internal calibration 
error to be measured. Fig.|Hlshows these repeat observa- 
tions for each run. Instead of plotting the difference be- 
tween the calibrated and standard index values for each 
star as above, in this figure the difference between cali- 
brated individual observations and the average of all ob- 
servations of that star is shown. 

Though not often listed, it is straightforward to calcu- 
late the uncertainty within the L ick/IDS system in any 
given observing run. I Wort he v et alJ \ 19941) lists the typi- 
cal error for an individual IDS measurement of each Lick 
index. Since each standard star was given equal weight 
in the calibration to the Lick/IDS system, the IDS in- 
ternal error was calculated as the average of the errors 
in the individual stars, weighted by the number of IDS 
observations. 
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The calibration to the Lick/IDS system is shown in 
Fig. The calibrated standard star observations from 
this work are in good agreement with both the Lick/IDS 
system and standard star observations from G93. 

2.7. Additional Errors 

Several sources of error not included in the simulations 
can also affect the equivalent width measurements to 
varying degrees. These include incorrect radial velocities, 
velocity dispersions, emission correction, and calibration 
error, as well as possible problems with the adopted data 
reduction procedures. 

Preliminary data reduction demonstrated that pub- 
lished radial velocities for the galaxies or standard stars 
do not necessarily agree with the observed spectra. The 
mismatches could be over 100 km/s in magnitude, with 
correspondingly large effects on the equivalent width 
measurements. This source of systematic error is avoided 
through the use of our empirical radial velocity measure- 
ments, as described above. Though the preliminary re- 
ductions did not show any similar problem with pub- 
lished velocity dispersions, empirical measurements were 
again chosen for the same reason. 

Velocity disper sion corrections are known to be very 
index-dependent ijTraeer et alJ ll998). The fitting func- 
tions from that work were used to translate a conser- 
vative error estimate a a ~ 10 km/s into index errors. 
These errors proved to be negligible for H/3 and Mgb. 
The actual errors in a are estimated to be ~ 3-5 km/s 
by the Pixfit code. The uncertainty in the measured ve- 
locity dispersion is the quadrature sum of the internal 
error measured by Pixfit with an estimate of the addi- 
tional uncertainty due to template mismatch. The latter 
quantity was estimated to be ~ 3 km/s based on compar- 
ing the measured velocity dispersions using several model 
spectra which fit the data nearly as well as the best fit 
spectrum. The measured velocity dispersions can eas- 
ily be compared with data in the literature, as shown in 
Fig- El The mean offset is 3 km/s with an RMS scatter 
of 11 km/s. 

The emission correction for H/3, H/y, and Fe5015 can 
drastically increase the error for those galaxies for which 
the correction is large. The error in [OIII]A5007 emis- 
sion was estimated by measuring the emission line in 
the difference spectrum (see above) using several differ- 
ent continuum levels. The average difference between 
these alternate continuum level equivalent width mea- 
surements and the zero continuum equivalent width mea- 
surement (which assumes a perfect match was achieved 
between the continuum levels of the model and ob- 
served spectra) was taken to be the uncertainty in the 
[OIIIJA5007 measurement. For many galaxies with little 
or no [OIII]A5007 emission, this uncertainty is negligi- 
ble, while for the galaxies with the largest [OIIIJA5007 
emission lines the uncertainty is so large as to render the 
emission corrections extremely suspect. The five galaxies 
with [OIIIJA5007 emission greater than 0.5A equivalent 
width have been omitted from subsequent analysis of the 
volume-limited sample due to the large uncertainties. 

An important potential source of systematic error is 
the method of background subtraction. Preliminary re- 
ductions tested two different techniques. The method 
using sky exposures is described above. A simpler tech- 
nique of simply subtracting a linear interpolation of the 



background between two user-defined regions near the 
edge of the slit was also tested. The latter technique is 
the mos t common in the li terature (G93 long slit obser- 
vations, [Beuinj^^ilj2002) , though the possibility of age 
and abundance gradients makes it dubious for galaxies 
of large angular size. For most galaxies in this study, 
the slit edges provide adequate sky subtraction. Reduc- 
ing the same galaxy using both sky subtraction meth- 
ods results in final index measurements which differ by 
much less than the total uncertainties. This check also 
indicates that the technique described above using the 
sky exposures does not introduce significant wavelength- 
dependent biases into the spectrum at the blue end. 

In addition to sky continuum subtraction, the bright 
Hg4358 night sky line contaminates the G4300 and H7F 
indices. The uncertainties of those indices are undoubt- 
edly larger than estimated, and the uncertainties are 
likely to vary with galaxy redshift depending on whether 
the sky line falls in the index bandpass or a sideband. 

If the input images have non-negligible subpixel offsets 
from one another in the dispersion direction, it is pos- 
sible that the effective instrumental resolution could be 
artificially degraded when the images are combined into 
the final two dimensional galaxy spectrum. To test for 
this, one galaxy (NGC 1209) was re-reduced in a dif- 
ferent manner following sky subtraction. Each input 
frame was reduced separately. Wavelength calibration 
was performed on the two dimensional images using stan- 
dard techniques (identify, reidentify, fitcoords, 
and transform in IRAF's longslit package). The 
r e /8 aperture on each was traced and extracted, then 
doppler corrected and smoothed to Lick/IDS resolution. 
The average index values showed no significant differ- 
ence from those derived using the previous data reduc- 
tion technique. Pixfit measured velocity dispersions less 
than 5 km/s smaller than those measured on the com- 
bined extraction previously. We conclude that there is 
little loss of effective instrumental resolution from com- 
bining input images. 

Since NGC 1209 was also used in the error simulations, 
the alternate reduction procedure above also provides an 
independent test of the estimated errors for each index. 
The alternate technique should measure only the random 
poisson errors between the spectra. Thus a priori one ex- 
pects smaller error bars than those estimated from the 
simulations which include several possible sources of sys- 
tematic error as well. The error comparison confirmed 
this expectation, except for Mgfe where the two error es- 
timates were equal. 

The final error bars adopted for the primary indices 
for each galaxy are a quadrature sum of the errors from 
the error simulation appropriate for that galaxy's S/N, 
the errors in a-correction and emission correction, and 
the internal and external calibration errors. Table^llists 
typical values for each of these quantities, as well as 
the uncertainty within the IDS system for a typical ob- 
serving run. Errors for the secondary indices are taken 
from the scatter in the measurements of each index in 
the NGC 1209 individual frames, added in quadrature 
as above with the a-correction, emission correction, and 
IDS calibration errors. Since the calibration error is a 
large contribution to many of these index errors, the re- 
sults should be reasonably accurate even for galaxies of 
significantly lower S/N than NGC 1209. 
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As a final consistency check, index measurements 
from this study are compared with measurements of the 
same galaxies from G93 and iDenicolo et al.l l)2005|) in 
Fig.HU An equivalent comparison between G93 and Deni- 
colo et al. is shown in Fig. 1 101 The measurements from 
this study are found to be in good agreement with G93. 
The Denicolo et al. sample is found to have larger error 
bars and systematic offsets in each measured quantity 
compared to both G93 and this study. Although crucial 
as an intermediate step in establishing the consistency 
between the latter two data sets, the Denicolo et al. sam- 
ple will not be included in the discussion in § 3. 

3. ANALYSIS 

3.1. Stellar Population Parameters 

The Balmer indices (H/3, ITyp, and H<5p) arc the 
most commonly used age indicators in stellar popula- 
tion analysis. The data in this study do not include the 
H^f index, but both H/ 3 and H7F have been measured. 
iKuntschner et ail l)2002fl presents a detail ed discussion of 
the re lative merits of each index (see also iCaldwell et aTl 
(2003) for a discussion of the H/3 index compared to 
higher order Balmer line ratios). Briefly, H/3 suffers from 
nearly twice the contamination of nebular emission lines 
(if present) filling in the stellar absorption. H7F also 
offers increased dynamic range between young and old 
stellar population models. Offsetting this, H/3 measure- 
ments normally achieve higher S/N due to the combi- 
nation of the red colors of early-type galaxies and the 
typical sensitivity curves of CCDs and spectrograph op- 
tics diminishing at bluer wavelengths. Both indices can 
suffer from m etal line contamination i n the continuum 
measurement. iStrader fe Brodie 1 l)2004|) showed that H/3 
suffers from additional variance caused by weak Fe I lines 
in the wavelength range of the index. As a result, ages 
based on H/3 can be less reliable than those based on H7 
or H5, as seen in the ext ragalactic globular c luster study 
of lPuzial l|200l (see also lHowell etall l|20?)l ). However, 
H7F can suffer from CH contamination and as noted in 
§ 2.7 is often affected by a bright night sky line, neither of 
which are issues for H/3. Pragmatically, H/3 is the age in- 
dicator of choice for this study, as the G93 data set does 
not include measurements of H7 or US. A comparison 
of age measurements using H/3 and H7F for the galaxies 
from this study will be presented below (§ 3.2). 

The fully calibrated and corrected equivalent widths 
for the galaxy sample are presented in Fig. ^] plotted 
on the H/3-[MgFe]' plane. The age-metallicity grid of 
the TMB models is also plotted. The composite index 
[MgFe]' was defined in TMB as an overall (i.e. [Z/H]) 
metallicity index that should be independent of [a/Fe] 
based on their calcula t ions using the response func- 
tions of lTripicco fc Belli l)1995|) . Members of the volume- 
limited sample observed by G93 are also plotted. The 
SSP age and metallicity for each galaxy was measured 
by linear interpolation among the TMB models in each 
quantity. 

The a-enhancement for each galaxy was measured us- 
ing the TMB models. Since age does not strongly affect 
the location of lines of constant [a/Fe] in the Mg&-(Fe) 
plane, each galaxy need not be compared individually to 
the model of precisely matching age as measured from 
Fig- [HI The 3 Gyr model was used to measure [a/Fe] 
for each galaxy, and an age correction was then applied. 



This age correction was calculated by linearly interpo- 
lating the offsets perpendicular to isoenhancement lines 
between models with ages of 3, 8, and 15 Gyr. The 
galaxy sample and the models are plotted in Fig. IT2"1 

Errors in age and metallicity are derived from the in- 
dex errors in H/3 and [MgFe]'. The extrema of the error 
ellipse defined by oh/3 and crj MgFe ]' in the directions per- 
pendicular to isochrone and isomctallicity lines were used 
to define at and az respectively. These stellar population 
uncertainties are highly correlated. Errors introduced 
by the interpolation of the TMB model are estimated 
to be 0.05-0.1 Gyr (depending on age) and < 0.01 dex 
in [Z/H], in each case negligible compared to the errors 
derived from the index uncertainties. 

The error in a-enhancement is derived in a similar fash- 
ion, with the added complication of including cr t . The 
extrema of the a^g^—a (Fc) error ellipse in the direction 
perpendicular to the isoenhancement lines are combined 
with ages of t ± cr t in the sense that maximizes the re- 
sulting spread in measured [a/Fe]. 

Parameter values and associated errors are presented 
in Table 0] Galaxies from G93 are listed in the bottom 
portion of the table. For consistency, the index measure- 
ments from G93 were used to determine stellar popula- 
tion para meters from TMB rath er than using the param- 
eters from lTrager et all l|2000aj) . 

3.2. Model Uncertainties 

The choice of model can strongly affect the values of 
the physical quantities one derives from spectral line in- 
dex measurements. The importance of measuring the 
G93 galaxies with respect to the TMB models is clearly 
shown in Figure 1131 The TMB models result in higher 
[a/Fe] (by an average of 0.06 dex) and [Z/H] (by an aver- 
age of 0.1 dex) than the Worthey models used by TFWG. 
This suggests that the earlier models yield the correct Fe 
abundance, but underestimate the a-elements and thus 
also [Z/H]. 

The H7 index is measured in the data from this work, 
and a comparison between the age estimates using H/3 
and H7 can be made as a check on the extent to which the 
choice of Balmer index affects the resulting age estimates. 
Within this Balmer line comparison, it is also instruc- 
tive to compa re the effect of us ing different SSP models. 
The models of lSchiavon I pfiol are used for this purpose. 
First, the Balmer index measurements are compared with 
one another in Fig. ^] With the exception of four out- 
liers, the two indices are tightly correlated as expected. 
FigurcslT31fc[THlshow the galaxies presented in this study 
plotted on the H/3 vs. (Fe) and H7F vs. (Fe) planes re- 
spectively. The TMB and S chiavon model grids are plot- 
ted on the former, and the lThomas. Maraston. fc KorrJ 
( 2004) and Schiavon model grids are plotted on the lat- 
ter. All models have [a/Fe] = 0. Ages derived us- 
ing each set of models are plotted against each other 
in Fig. 1171 As is apparent from a visual inspection of 
the grids, an age offset of several Gyr exists between the 
Thomas et al. and Schiavon model grids, in the sense of 
the Thomas et al. models yielding (on average) younger 
ages on H/3 grids and older ages on H7F grids compared 
to the Schiavon models. The outliers in Fig. 1141 remain 
outliers in Fig. El using either set of models. The quan- 
titative difference between age estimates from the two 
sets of models suggests that the choice of SSP model is 
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more important than the choice of Balmer index when 
attempting to measure ages by this method. Since the 
two sets of models each produce consistent age measure- 
ments for some range of ages (Thomas et al. models are 
consistent for old galaxies; Schiavon models are consis- 
tent for young galaxies), for the purposes of this study 
there is no clear reason to prefer one to the other. There- 
fore the choice to use the TMB models and measure SSP 
parameters on model grids with H/3 as the age-sensitive 
index does not introduce any additional systematic error 
due to choice of models or Balmer index. Note that the 
above comparisons do indicate the presence of system- 
atic errors arising from both sources, however the same 
would be true regardless of the models or Balmer index 
selected. Also, an inconsistency between ages derived us- 
ing H/3 and H7 does not necessarily imply that the SSP 
model being used is incorrect in some way. As noted in 
§ 1, single stellar populations are a necessary simplifying 
assumption, but galaxies in general will have stellar pop- 
ulations of varying ages and metallicities. Such a com- 
posite stellar population can easily yield different ages 
derived using H7 as the age sensitive index compared to 
using H/3 as the age sensitive index. 

3.3. Multiple Stellar Population Analysis 

Current stellar population analysis techniques rely on 
the single stellar population (SSP) assumption, that a 
single age, metallicity, and a-element abundance ratio 
can characterize the stellar population of a galaxy. More 
realistic stellar population models can be constructed, 
however there is no widely applicable method to fully 
constrain the problem — even in the simplest case of as- 
suming two starbursts instead of one, the number of free 
parameters is more than double that of the SSP case 
(two ages, two [Z/H], two [a/Fe], and the mass ratio 
between the two bursts). Allowing for a larger number 
of bursts or continuous star formation further undercon- 
strains the problem. In a few galaxies, SSP observations 
of their globular cluster (GC) populations (for which the 
SSP assumption is correct rather than a means to ob- 
tain a luminosity-weighted average) may allow two-burst 
multiple stellar population models. Almost all ellipti- 
cal galaxies contain at least two populat ions of GCs, as 
seen in bimodal color distributions (e. g. lZepf fc Ashmanl 
1993). Measurements of these two GC populations al- 
low the two sets of stellar population parameters for a 
two-burst ga l axy m odel to be set in a consistent way. 
iHowell et alJ (|2004T ) showed that this does not always 
work; if the youngest GC population is coeval with the 
SSP age of the galaxy the resulting mass ratio between 
the two starbursts will necessarily be strongly weighted 
toward s the younger burst. As shown in iTrager et al.l 
( 2000b) , for example, although a small mass fraction of 
young stars is sufficient to change the measured SSP age 
by several Gyr, mass fractions of 50% or more are re- 
quired to produce SSP ages consistent with that of the 
young population itself. Such large burst fractions are 
contrary to most expectations for galaxy formation. 

The GC populations of a few galaxies in the volume- 
limited sample have been studied spectroscopically. 
With stellar population information for two GC popu- 
lations, the age, [Z/H], and [a/Fe] for each burst can be 
estimated rather than being free parameters. In prac- 
tice, [a/Fe] can be ignored as the best available models 



for multiple ste llar population analysis remain those of 
iWorthevI (|1QQ4T ) via Worthey's Dial-a-Galaxy web page 3 , 
which do not account for non-solar abundance ratios. 

The requisite data exist for four galaxies: NGC 3115, 
NGC 3610, NGC 4365, a nd NGC 584 6 NG C 3610 
is discussed in de tail in IHowell et alJ l)2004(K Us- 
ing IR photometry, iPuzia et al.l l)2002fl found candidate 
intermediate-age GCs associated with NG C 4365. These 
candid ates were seemingly confirmed by lLarsen et alJ 
(2003) finding a young and metal rich GC population 
with an age between 2 -5 Gyr and approxi mately solar 
metallicity. However, iBrodie et alJ l)2005t) found that 
based on all Balmer lines the GC population was coeval 
within the errors, w ith an age > 10 Gyr. The data from 
lLarsen et al.l l)2003() were biased by low S/N in the sense 
that some but no t all of the intermediate age candidates 
from IPuzia et alJ l)2002j) appeared young using only the 
H/3 index. Since the GC population is older than the 
SSP age of the galaxy, no viable two burst model can 
be constructed using the GC populations to constrain 
the burs ts. NGC 3115 appears t o be a similar case, with 
old GCs llKuntsch ner et a l.l2002D around an intermediate 
age galaxy. ITrager et al.l l)1998fi suggest that the galaxy 
itself is also old, however the NGC 3115 data presented 
here are of exemplary S/N with negligible emission cor- 
rection, and conclusively indicate an age between 4- 
7 Gyr depending on whether H/3 or H7 is used as the age- 
sensitive index. The aperture used in this work i s large r 
tha n that of the data analyz ed in ITrager et alJ |l998), 
and lBassin fc Bonattol l)2003l) showed that the contribu- 
tion of a young stellar population increases with radius 
on the scale covered by aperture used here. The discrep- 
a ncy between the mea surement presented here and that 
of ITrager et al.1 i|1998fl is therefore explained by differing 
stellar pop ulations in the ap ertures used in each study. 
NGC 5846 l|Puzia et all2004ft has uniformly old GC pop- 
ulations at a range of metallicities. Since the galaxy has 
an old SSP age (Table 0J), the SSP assumption appears 
to be a reasonable approximation to the star formation 
histories of this galaxy. 

3.4. Correlations Between SSP Parameters 

In Fig. 1181 all galaxies in the volume-limited sample 
observed thus far are plotted on the metallicity "hyper- 
plane" of TFWG. Error ellipses in the hyper-plane pro- 
jections were calculated from the errors in the individual 
indices. Index errors were treated as independent quan- 
tities, as sky subtraction error is negligible at the centers 
of galaxies. Monte Carlo simulations were run propagat- 
ing index error realizations for a particular galaxy into 
physical parameters and thence to the principal compo- 
nent axes of the hyper-plane. The RMS scatter about 
the major and minor axes of the Monte Carlo distri- 
bution in the hyper-plane was used for the major and 
minor axis lengths of the la error ellipses. The actual 
error distribution is not quite symmetrical: as one ex- 
pects from looking at the model grids, a fixed level of 
index scatter leads to unequal parameter scatter in the 
young and metal-rich directions compared to the old and 
metal-poor directions. The hyper-plane is clearly due to 
physical variations between galaxies, not measurement 
error or bias. The volume-limited and G93 samples have 
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similar distributions within the hyper-plane. 

The Z-plane from TFWG (Fig.[jj"J) provides an intrigu- 
ing challenge to theories of galaxy formation. This pro- 
jection shows that a linear relationship exists between 
[Z/H], log cr, and log t. TFWG showed that no other 
galaxy parameters significantly affected this fit. The fit 
from TFWG is plotted; however a slight offset would 
better fit the data presented here. This is the result of 
using the TMB SSP models instead of the Worthey mod- 
els used by TFWG; recall that Fig. If 31 showed a constant 
offset between the two sets of models in the [Z/H] mea- 
surements. This offset in [Z/H] requires a corresponding 
offset in the Z-plane fit. As discussed in detail in TFWG, 
single burst models cannot plausibly produce a Z-plane 
relation — a linear relation at the present epoch would 
be noticeably curved at past and future times, as old 
objects change less in log t than young objects. TFWG 
showed that a two (or several) burst "frosting" model can 
in principle maintain the Z-plane over time. Such frost- 
ing populations most plausibly form from gas within the 
galaxy and enriched by previous star-formation events; it 
is unlikely that accreted gas or a merger event would have 
the appropriate metallicity to keep the galaxy on the Z- 
plane. However, the young and metal-rich end of the Z- 
plane tends to be populated by likely merger remnants 
such as NGC 36 f 0, the most extreme such galaxy in the 
volume-limited sample. As shown in Fig. [""""| NGC 36f 
lies directly on the Z-plane, as do other similar galax- 
ies from the G93 sample which are also believed to be 
merger remnants (TFWG). 

The volume-limited sample is compared with the G93 
sample in all four principal component axes, as well as 
position along the edge-on projection of the Z-plane. K-S 
tests indicate that the volume-limited sample is consis- 
tent with being drawn from the same distribution as the 
G93 sample in all five dimensions. This confirms that the 
G93 galaxies are representative of the early-type galaxy 
population in the local universe. 

3.5. Correlations with Other Parameters 

Having derived stellar population parameters for the 
galaxies in the volume-limited sample, as well as rederiv- 
ing the stellar population parameters for the G93 sam- 
ple using the same models, we now turn our attention 
to other physical parameters of these galaxies. Figure l""0l 
shows the stellar population parameters age, [Z/H], and 
[a/Fe] plotted against three structural parameters, cr, 
Mr, and ry,, as well as the anisotropy parameter (v /a) * 
'Bender. Burstein. fc Fabellf 992tlFaber ^et alll997*) . the 
core profile slope "y llLauer et al.1 2005), and the B — 
V co lor gradient ijldiart. Mi chard & de Freitas P achecol 
2002). These quantities are listed in Table "5J Note that 
the relatively small magnitude range covered by this sam- 
ple limits the ability to investigate correlations involv- 
ing luminosity or quantities closely tied to luminosity. 
For each pair of parameters plotted against each other, 
the correlation statistic was calculated, excluding those 
galaxies with large emission corrections. The Kendall's 
t rank correlation test was used to identify correlations 
which were not significant. Unfortunately the Kendall's 
r test can only rule out correlations; it cannot iden- 
tify which correlations are significant, as demonstrated 
in trials where it found nominally significant correlations 
between NGC number and various physical parameters. 



The Spearman rank correlation test was found to have 
similar behavior. Thus a correlation is judged to be sig- 
nificant if it is not ruled out by the Kendall's r test and 
the absolute value of the correlation statistic is greater 
than a threshold value, arbitrarily chosen to be 0.3. By 
these criteria, significant correlations are found between 
a and both age and [a/Fe]; and between log r e and age. 
The relation between [a/Fe] and a is a representation of 
the well-known Mg-cr relation. The other two relations 
indicate that the smallest galaxies, measured either by 
radius or mass, are younger than larger galaxies, as ex- 
pe cted from the "dow nsizing" galaxy formation scheme 
of|C owie et all (GUI. This result is in agreement with 
the more c omprehensive study of low velocity dispersion 
galaxies bv lCaldwell et alJ l|2003|) . 

Early-type galaxies lie on planes in two different pa- 
rameter spaces, the Fundamental Plane and the metal- 
licity hyper-plane. The existence of a mapping from one 
plane to the other would imply a relation connecting the 
size and brightness of a galaxy to its age and chemi- 
cal composition. To search for such a mapping, the re- 
space parametCTi^tion of the Fundamental Plane was 
used (Ben der et alJlT992|) . The two planes are shown in 
Fig. EI] The IRAF task GEOMAP found no combination 
of scale factors, rotations, and translations that can map 
the K\-K2 face-on projection of the Fundamental Plane 
into the PC1-PC2 face-on projection of the metallicity 
hyperplane. The best fit resulted in average offset dis- 
tances of 0.15 in the PC1-PC2 plane (average offset of 
0.07 in the PCI axis and 0.11 in the PC2 axis individu- 
ally). The largest offsets occured for the most extreme 
galaxies in one or both axes, up to a difference of 0.52 in 
PC2 for one galaxy. As a result, the transformed Funda- 
mental Plane distribution covers a smaller portion of the 
PC1-PC2 plane than the metallicity hyper-plane. The 
transformed Fundamental Plane distribution is found to 
be consistent with the metallicity hyper-plane by K-S 
tests, though the consistency is marginal (K-S probabil- 
ity of 8%) on the PC2 axis. Also, as previously noted, 
the sample presented here covers a relatively small range 
of luminosity and therefore may lack leverage to reveal 
important details within these parameter spaces. 

The stellar brightness profiles of 27 of the galaxies in 
both the volume- limited s ample and the G93 sample have 
been studied with HST (jFaber et alJll997t lLauer et aT] 
2005). Figure 12*01 showed no significant correlations be- 
tween the brightness profile slope 7 and any stellar pop- 
ulation parameter, though this test was restricted to 
the subset of the volume-limited sample for which HST 
surface brightness profiles are available. Combining all 
available galaxies, including G93 galaxies not included 
in the volume-limited sample, the distribution of cr, age, 
[Z/H], and [a/Fe] among both core and power-law galax- 
ies was examined using K-S tests. The two classes of 
galaxies are consistent with having the same distribution 
of age and metallicity (probabilities 0.087 and 0.82, re- 
spectively), while core galaxies have significantly higher 
velocity dispersions and [a/Fe] values than power-law 
galaxies (probabilities of 0.019 and 0. 036, respectively) . 
This is as expected given the result of (jFaber et al.lfl 9971 
that the most luminous galaxies all have core profiles, 
while all galaxies below a certain luminosity threshold 
(My > —20.5) have power-law profiles. To further ex- 
plore the meaning of these relations, the distribution 
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of stellar population parameters between g alaxies with 
boxy o r disky isophotes was examined. Since lFaber et al.1 
( 1997) found a strong correlation between galaxies with 
disky isophotes and power law core profiles, the expec- 
tation is that the disky galaxies will have lower [a/Fe] 
and a measurements. I sophotal shape me asurements for 
15 galaxies are li sted inlldiart et all (l2002h . and the 04/0 
parameter from iFaber et"ldTl(T99^ was used to deter- 
mine the shapes of an additional 13 galaxies. Boxy and 
disky galaxies are consistent with being drawn from the 
same distribution in cr, age, and [Z/H] (probabilities of 
0.43, 0.23, and 0.75, respectively), while disky systems 
have significantly smaller [a/Fe] than boxy systems as 

e xpected (probability of 0. 048). 

iCarollo et all <|1997j) and lSpaans fc Carollol <|1997f) ar- 
gue that more massive, spherical systems form stars 
earlier and for a shorter duration t han less massive, 
more fl attened systems. The models of lSoaans fc Carollol 
(1997) show that flattened systems lose metals more eas- 
ily due to supernova blow-out, which has the consequence 
of delaying metal enrichment in the outer regions and 
thus delaying the peak of the star formation rate. This 
peak occurs when a metal-enriched multiphase ISM has 
built up, at which time gas can cool (and thus form stars) 
much more efficiently. By this argument, flattened galax- 
ies — disky, power law profile, high (v/a)*, low mass — 
should have lower [a/Fe] and younger ages than boxy, 
core profile, low (w/u)», high mass galaxies. Many of 
these predictions are confirmed, as mentioned above: a 
significant correlation is seen in the sense that smaller, 
less massive galaxies are younger than larger, more mas- 
sive galaxies fFig. I20|) ; galaxies with power-law profiles 
have significantly lower [a/Fe] than galaxies with core 
profiles; and galaxies with disky isophotes have signifi- 
cantly lower [a/Fe] than galaxies with boxy isophotes. 
No statistically significant difference in age distribution 
is seen between core and power-law galaxies or boxy and 
disky galaxies, however, nor are any correlations between 
stellar population parameters and (u/er)* found to be sig- 
nificant. 

The models of iSoaans fc CaroTIol l)1997ft t reat galaxies 
as evo lving and forming stars in isolation. As lFaber et al.l 
( 1997) summarized, gas-rich merger events lead natu- 
rally to disk y galaxies with power-law profiles. Simu- 
lations (e. g. iBarnes fc Hernauistl fl996|> show that gas 
can destroy box orbits and that the gas (and thus star 
formation) in th e remnant object is centrally concen- 
trated. However. IFaber et alJ l)1997[) also noted that the 
merger of the supermassive black holes from the pro- 
genitor galaxies can produce core brightness profiles by 
ejecting stars from the center of the merger remnant. It 
is therefore unclear to what extent major mergers are ex- 
pected to yield correlations between structural and stel- 
lar population parameters. 

To search for additional correlations, principal compo- 
nent analyses were run on several combinations of pa- 
rameters (Table EJ. The input parameters were each re- 
duced to zero mean and unit variance prior to the anal- 
ysis. Combinations of parameters including luminosity 
(My), one structural variable (7 or the break radius r&), 
and one SSP variable were studied; in addition, all seven 
galaxy paramet ers were ana lyzed together. Using the 
Kaiser criterion (Kaiser 1960), only those principal com- 
ponents which account for at least as much variance as 



a single normalized input variable are listed in Table 
No new planes were discovered in which two orthogonal 
linear combinations of the input variables could explain 
almost all of the variance in the three dimensional input 
parameter space. The distributions of galaxies in the pa- 
rameter spaces studied by this method are triaxial, in 
some cases close to prolate spheroids and in other cases 
close to oblate spheroids. 



4. SUMMARY AND CONCLUSIONS 

A new set of spectra of nearby elliptical galaxies has 
been obtained. Combining this new data set with a se- 
lected portion of the G93 sample, a volume-limited sam- 
ple of 45 galaxies has been defined, of which 35 have been 
observed. The ten missing galaxies have a distribution 
in luminosity and velocity dispersion that is similar to 
the full set of 45 galaxies. 

It is worth noting that there is no significant qualita- 
tive difference between SSP models; for almost all galax- 
ies the use of different models produces only a systematic 
offset in stellar population parameters. Although differ- 
ent choices of models will result in different absolute SSP 
quantities, comparing one galaxy to another within the 
same set of models yields consistent results independent 
of the choice of models. This quantitative offset between 
different models is greatest when higher order Balmer 
lines such as H7 are used as the age-sensitive index in 
age-metallicity grids. The difference in ages derived us- 
ing different sets of models is found to be greater than the 
difference in ages derived using different Balmer indices 
in model grids. 

A method to constrain the input parameters for a two- 
burst multiple stellar population analysis has been pro- 
posed. However, stellar population data are available for 
the GC systems of only a few galaxies; four such galaxies 
are investigated here. In none of these cases is a two- 
burst model based on GC age measurements viable. Two 
galaxies are coeval with their GC populations, suggest- 
ing that SSP analysis is adequate for those galaxies. The 
other two galaxies have SSP ages significantly younger 
than any of the associated GCs. 

The existence of the metallicity hyper-plane and the 
Z-plane of TFWG is confirmed by this new data set. 
The distribution of galaxies within these planes is very 
similar to that shown by TFWG, and this distribution is 
shown to be due to physical variations well in excess of 
measurement uncertainties. 

The G93 galaxy sample analyzed by TFWG is repre- 
sentative of the local early-type galaxy population. This 
work is in good statistical agreement with their results. 

Comparisons between stellar population parameters 
and indicators of a g alaxy's internal structure support 
the idea introduced bv|S paans fc Carollol l|1997fl that size 
and angular momentum (i.e. large and spherical com- 
pared to smaller and flattened) are key distinctions be- 
tween power law and core profile galaxies. Disky, power- 
law profile galaxies are found to have significantly lower 
[a/Fe] than boxy, core profile galaxies, though no signif- 
icant trend in age is seen between core and power-law or 
boxy and disky galaxies. Age is found to decrease with 
size, measured either by radius or mass, as expected in 
downsizing models for galaxy formation. 
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TABLE 1 
Volume Limited Galaxy Sample 



Name 


a(J2000) 


<5(J2000) 


S/N 


Observed 


NGC 


584 


01 h 31 m 20.7 s 


-06°52'06" 




G93 


NGC 


596 


01 h 32 m 52.08 s 


-07°01'54.6" 


147.0 


2002 November 3 


NGC 


720 


01 h 53 m 00.4 s 


_13°44'18" 




G93 


NGC 


821 


02 h 08 m 21.0 s 


+10°59'44" 




G93 


NGC 


1052 


02 h 41 m 04.80 B 


-08°15'20.8" 


205.8 


2000 November 26 


NGC 


1172 


03 h 01 m 36.0 s 


-14°50'11" 


97.0 


2002 November 3 


NGC 


1199 


03 h 03 m 38.6 s 


-15°36'51" 




Not Observed 


NGC 


1209 


03 h 06 m 03.0 s 


-15°36'40" 


149.2 


2000 November 27 


NGC 


1400 


03 h 39 m 31.0 s 


-18°41'22" 


121.0 


2003 January 2 


NGC 


1407 


03 h 40 m 11.8 s 


-18° 34' 48" 


164.3 


2000 November 26 


NGC 


2768 


09 h ll m 37.50 s 


+60°02'15.0" 


89.7 


2002 November 4 


NGC 


2974 


09 h 42 m 32.96 s 


-03°41'55.2" 




Not Observed 


NGC 


3115 


10 h 05 m 13.42 s 


-07°43'06.5" 


305.0 


2001 March 28 


NGC 


3156 


10 h 12 m 41.08 s 


+03°07'50.2" 




Not Observed 


NGC 


3193 


10 h 18 m 24.88 s 


+21°53'38.6" 


108.9 


2004 March 27 


NGC 


3226 


10 h 23 m 27.00 s 


+19°53'54.4" 


101.9 


2001 March 25 


NGC 


3607 


ll h 16 m 54.08 s 


+18°03'11.6" 


141.0 


2003 January 2 


NGC 


3608 


ll h 16 m 58.7 s 


+18°08'57" 




G93 


NGC 


3610 


ll h 18 m 25.83 s 


+58°47'13.6" 


136.4 


2001 March 27 


NGC 


3613 


ll h 18 m 36.12 s 


+58°00'04.5" 


107.4 


2001 March 27 


NGC 


3640 


ll h 21 m 06.74 s 


+03° 14' 08.1" 


151.9 


2001 March 28 


NGC 


3962 


ll h 54 m 40.0 s 


-13°58'30" 




Not Observed 


NGC 


4125 


12 h 08 m 05.71 s 


+65° 10' 24.5" 




Not Observed 


NGC 


4168 


12 h 12 m 16.9 s 


+13°12'20" 


103.7 


2001 March 25 


NGC 


4278 


12 h 20 m 06.82 s 


+29°16'50.7" 




Not Observed 


NGC 


4365 


12 h 24 m 27.87 s 


+07°19'04.9" 


158.3 


2001 March 26 


NGC 


4374 


12 h 25 m 03.7 s 


+12°53'14" 




G93 


NGC 


4406 


12 h 26 m 11.74 s 


+12°56'46.4" 




Not Observed 


NGC 


4472 


12 h 29 m 46.5 s 


+07°59'48" 




G93 


NGC 


4473 


12 h 29 m 48.87 s 


+13°25'45.7" 


186.7 


2001 March 28 


NGC 


4486 


12 h 30 m 49.42 s 


+12°23'28.0" 


150.9 


2004 March 27 


NGC 


4552 


12 h 35 m 39.9 s 


+12°33'55" 




G93 


NGC 


4621 


12 h 42 m 02.49 s 


+11°38'48.7" 


166.3 


2004 March 27 


NGC 


4636 


12 h 42 m 49.7 s 


+02°41'18.4" 




Not Observed 


NGC 


4649 


12 h 43 m 39.7 s 


+11°33'09" 




G93 


NGC 


4697 


12 h 48 m 35.8 s 


-05° 48' 00" 




G93 


NGC 


5322 


13 h 49 m 15.19 s 


+60°11'26.2" 




Not Observed 


NGC 


5485 


14 h 07 m 11.5 s 


+55°00'07" 




Not Observed 


NGC 


5576 


14 h 21 m 04.11 s 


+03°16'13.5" 


144.0 


2001 March 26 


NGC 


5638 


14 h 29 m 40.4 s 


+03° 14' 04" 




G93 


NGC 


5812 


15 h 00 m 57.0 s 


-07°27'19" 




G93 


NGC 


5813 


15 h 01 m 11.2 s 


+01°42'08" 




G93 


NGC 


5831 


15 h 04 m 07.2 s 


+01°13'15" 




G93 


NGC 


5846 


15 h 06 m 29.3 s 


+01°36'21" 


153.2 


G93, 2001 March 28 


NGC 


6703 


18 h 47 m 18.9 s 


+45°33'02" 




G93 



TABLE 2 

Index Measurements: r e /8 Aperture 



Galaxy a Ca4227 G4300 H7 F Fc4383 Ca4455 Fc4531 C 2 4668 H/3 [OIII] Fc5015 Mg 2 Mgb Fe5270 Fe5335 Fc5406 [MgFc] ' 



NGC 596 
NGC 1052 
NGC 1172 
NGC 1209 
NGC 1400 
NGC 1407 
NGC 2768 
NGC 3115 
NGC 3193 
NGC 3226 
NGC 3607 
NGC 3610 
NGC 3613 
NGC 3640 
NGC 4168 
NGC 4365 
NGC 4473 
NGC 4486 
NGC 4621 
NGC 5576 
NGC 5846 



151 
4 

215 
4 

113 

4 

225 
4 

285 
5 

296 
4 

211 

4 

276 
4 

228 
4 

180 
5 

231 
4 

172 
4 

221 
4 

178 
4 

179 
5 

270 
4 

201 
4 

371 
5 

260 
4 

190 
4 

243 
4 



1.19 
0.10 
1.20 
0.10 
1.21 
0.10 
1.39 
0.10 
1.39 
0.10 
1.47 
0.10 
1.46 
0.09 
1.55 
0.10 
1.27 
0.09 
1.09 
0.10 
1.14 
0.10 
0.89 
0.10 
1.29 
0.10 
1.20 
0.10 
1.08 
0.10 
1.58 
0.10 
1.31 
0.10 
1.71 
0.09 
1.45 
0.09 
1.10 
0.10 
1.40 
0.10 



5.33 
0.14 
5.72 
0.14 
5.29 
0.14 
5.35 
0.14 
5.08 
0.14 
4.80 
0.14 
5.36 
0.18 
5.57 
0.14 
5.17 
0.18 
5.85 
0.16 
5.27 
0.14 
4.99 
0.16 
5.16 
0.16 
5.22 
0.16 
5.07 
0.16 
5.41 
0.16 
5.62 
0.16 
5.47 
0.18 
5.50 
0.18 
5.36 
0.16 
5.09 
0.16 



-1.12 

0.13 

-2.17 

0.38 

-1.16 

0.15 

-1.85 

0.12 

-1.80 

0.14 

-1.18 

0.11 

-1.45 

0.20 

-1.78 

0.11 

-1.64 

0.13 

-2.17 

0.12 

-1.53 

0.14 

-0.52 

0.09 

-1.46 

0.08 

-1.29 

0.08 

-1.08 

0.08 

-1.48 

0.08 

-2.01 

0.08 

1.37 

0.16 

-2.07 

0.11 

-1.04 

0.08 

-1.16 

0.10 



4.54 
0.21 
6.50 
0.09 
4.38 
0.21 
5.58 
0.09 
4.72 
0.21 
6.29 
0.09 
5.36 
0.20 
5.61 
0.09 
5.10 
0.20 
4.75 
0.18 
4.86 
0.21 
4.72 
0.18 
5.00 
0.18 
4.88 
0.18 
4.85 
0.18 
5.24 
0.18 
5.48 
0.18 
5.71 
0.20 
5.34 
0.20 
4.66 
0.18 
5.39 
0.18 



2.06 
0.14 
1.43 
0.05 
1.35 
0.14 
1.58 
0.05 
1.47 
0.14 
1.30 
0.05 
1.50 
0.05 
1.53 
0.05 
1.59 
0.05 
1.39 
0.07 
1.47 
0.14 
1.52 
0.07 
1.77 
0.07 
1.51 
0.07 
1.59 
0.07 
1.95 
0.07 
1.70 
0.07 
1.30 
0.05 
1.54 
0.05 
1.41 
0.07 
1.51 
0.07 



3.62 
0.09 
3.78 
0.11 
3.36 
0.09 
3.76 
0.11 
3.90 
0.09 
3.87 
0.11 
3.80 
0.11 
4.03 
0.11 
3.94 
0.11 
3.67 
0.07 
3.78 
0.09 
3.57 
0.07 
3.35 
0.07 
3.61 
0.07 
3.52 
0.07 
4.04 
0.07 
3.79 
0.07 
4.24 
0.11 
4.03 
0.11 
3.57 
0.07 
3.70 
0.07 



6.60 
0.20 
8.24 
0.12 
4.97 
0.20 
7.75 
0.12 
7.47 
0.20 
8.80 
0.12 
5.23 
0.11 
8.52 
0.12 
7.52 
0.11 
7.26 
0.19 
8.22 
0.20 
8.08 
0.19 
7.43 
0.19 
6.83 
0.19 
6.65 
0.19 
8.78 
0.19 
7.64 
0.19 
8.71 
0.11 
8.39 
0.11 
7.32 
0.19 
8.12 
0.19 



1.94 

0.06 

1.21 

0.7 

1.94 

0.15 

1.27 

0.12 

1.33 

0.11 

1.40 

0.06 

1.62 

0.35 

1.66 

0.04 

1.50 

0.16 

1.73 

0.17 

1.53 

0.11 

2.39 

0.08 

1.85 

0.06 

1.87 

0.06 

1.87 

0.07 

1.55 

0.06 

1.74 

0.06 

1.07 

0.25 

1.26 

0.10 

2.10 

0.06 

1.44 

0.12 



-0.06 

0.03 

-3.71 

1.0 

-0.65 

0.20 

-0.23 

0.15 

-0.21 

0.13 

-0.01 

0.02 

-0.58 

0.47 

-0.03 

0.01 

-0.45 

0.18 

-1.48 

0.23 

-0.30 

0.13 

-0.03 

0.08 

0.00 

0.01 

-0.03 

0.02 

-0.26 

0.05 

-0.01 

0.01 

-0.01 

0.01 

-.57 

0.32 

0.0 

0.01 

-0.01 

0.02 

-0.23 

0.15 



5.43 
0.28 
1.90 
0.64 
4.91 
0.30 
5.91 
0.22 
6.03 
0.29 
5.90 
0.20 
4.29 
0.58 
6.23 
0.20 
5.28 
0.58 
4.14 
0.18 
5.89 
0.28 
5.51 
0.12 
5.92 
0.11 
5.79 
0.11 
5.06 
0.11 
5.57 
0.11 
5.73 
0.11 
3.63 
0.58 
5.66 
0.58 
5.52 
0.11 
4.81 
0.11 



0.260 
0.003 
0.340 
0.003 
0.238 
0.003 
0.326 
0.003 
0.336 
0.003 
0.350 
0.003 
0.200 
0.003 
0.325 
0.003 
0.318 
0.003 
0.294 
0.003 
0.313 
0.003 
0.255 
0.003 
0.289 
0.003 
0.266 
0.003 
0.258 
0.003 
0.330 
0.003 
0.315 
0.003 
0.348 
0.003 
0.348 
0.003 
0.271 
0.003 
0.333 
0.003 



3.89 
0.05 
5.96 
0.05 
3.89 
0.06 
4.99 
0.05 
5.32 
0.06 
5.41 
0.05 
4.36 
0.07 
5.04 
0.05 
4.61 
0.07 
4.73 
0.06 
4.76 
0.05 
4.04 
0.05 
4.47 
0.06 
3.96 
0.05 
3.98 
0.06 
5.18 
0.05 
4.73 
0.05 
5.78 
0.06 
5.09 
0.06 
4.13 
0.05 
5.20 
0.05 



2.91 
0.06 
3.05 
0.06 
2.70 
0.09 
3.32 
0.06 
3.01 
0.09 
3.50 
0.06 
2.81 
0.10 
3.33 
0.05 
3.02 
0.10 
2.78 
0.09 
2.96 
0.06 
2.77 
0.06 
3.21 
0.09 
2.96 
0.06 
2.92 
0.09 
3.32 
0.06 
3.22 
0.06 
2.88 
0.07 
3.15 
0.07 
3.03 
0.06 
3.13 
0.06 



2.46 
0.09 
2.78 
0.09 
2.27 
0.13 
2.78 
0.09 
2.75 
0.13 
2.83 
0.09 
2.70 
0.14 
3.21 
0.07 
2.67 
0.14 
2.77 
0.13 
2.78 
0.09 
3.10 
0.09 
3.15 
0.13 
2.84 
0.09 
2.69 
0.13 
3.19 
0.09 
3.02 
0.09 
3.22 
0.11 
3.07 
0.11 
2.83 
0.09 
3.03 
0.09 



1.72 
0.09 
1.88 
0.09 
1.61 
0.09 
1.72 
0.09 
1.82 
0.09 
1.86 
0.09 
1.92 
0.09 
2.07 
0.09 
2.16 
0.09 
1.78 
0.09 
1.83 
0.09 
1.81 
0.09 
2.00 
0.09 
1.80 
0.09 
1.61 
0.09 
2.07 
0.09 
1.81 
0.09 
2.37 
0.09 
2.15 
0.09 
1.79 
0.09 
1.91 
0.09 



3.29 
0.05 
4.21 
0.05 
3.17 
0.07 
3.98 
0.05 
3.95 
0.08 
4.23 
0.05 
3.48 
0.08 
4.08 
0.04 
3.67 
0.08 
3.62 
0.08 
3.72 
0.05 
3.40 
0.05 
3.78 
0.07 
3.40 
0.05 
3.37 
0.07 
4.12 
0.05 
3.87 
0.05 
4.15 
0.07 
3.99 
0.07 
3.50 
0.05 
4.02 
0.05 
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TABLE 3 
Primary Index Error Table 



Error 


H/3 


Mgb 


Fc5270 


Fe5335 


Simulations 


0.05 


0.03 


0.04 


0.06 


Velocity Dispersion 


0.00 


0.00 


0.01 


0.02 


Calibration - External 


0.05 


0.05 


0.04 


0.06 


Calibration - Internal 


0.02 


0.02 


0.02 


0.04 


Lick/IDS System 


0.07 


0.07 


0.08 


0.08 
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TABLE 4 
Stellar Population Parameters 



Galaxy 






Age 




[Z/H] 




[a/Fe] 


a a 




(km/s) 


(km/s) 


(Gyr) 


(Gyr) 


(dex) 


(dex) 


(dex) 


(dex) 






This work, 


TMB models: 








in uv.' uuyu 


1 ^1 


4 


4.4 


0.7 


0.22 


U. 1 


n 1 q 


u.ui 




215 


4 


16 


14 


0.42 


0.3 


0.44 


0.12 


wrp 1 i 79 

IN v_i v_/ 11 ( Z 




4 


4.8 


1.5 


0.13 


n iq 

U.lo 


u.zu 


u.uo 


in vjL' izuy 


225 


4 


15.6 


3.0 


0.28 


u.uo 


U.ZJ 


u.uo 


wnp 1400 


285 


5 


14.2 


3.0 


0.31 


0.13 


0.35 


0.04 


NOC 1 407 


296 


4 


9.5 


2.2 


0.56 


0.07 


0.30 


0.04 


IN v_ I V_ £ ( Uo 


211 


4 


10 


7.0 


0.14 


0.25 


0.22 


0.10 


IN VI V_/ OllU 


276 


4 


3.9 


0.7 


0.65 


0.06 


0.25 


0.03 


in oiyo 


228 


4 


11.8 


3.2 


0.20 


12 


24 


u.uu 


IN uv.' OZiZiU 


180 


5 


6.1 


3.5 


0.33 


0.17 


0.32 


0.06 


NGC 3607 


231 


4 


10.6 


2.3 


0.27 


0.1 


0.27 


0.04 


NGC 3610 


172 


4 


1.7 


0.1 


0.76 


0.16 


0.28 


0.03 


NGC 3613 


221 


4 


3.3 


1.0 


0.59 


0.12 


0.17 


0.05 


NGC 3640 


178 


4 


4.9 


1.0 


0.26 


0.08 


0.13 


0.04 


NGC 4168 


179 


5 


5.0 


1.4 


0.24 


0.1 


0.17 


0.05 


NGC 4365 


270 


4 


5.9 


1.9 


0.59 


0.08 


0.26 


0.04 


NGC 4473 


201 


4 


4.0 


1.3 


0.56 


0.08 


0.24 


0.04 


NGC 4486 


371 


5 


19.6 


7.5 


0.27 


0.18 


0.36 


0.08 


NGC 4621 


260 


4 


15.8 


4.0 


0.29 


0.11 


0.23 


0.05 


NGC 5576 


190 


4 


2.5 


0.4 


0.60 


0.11 


0.21 


0.05 


NGC 5846 


243 


4 


10.6 


4.4 


0.44 


0.14 


0.28 


0.03 



Galaxies from G93, volume-limited sample, TMB models: 



NGC 0584 


236 


3 


2.4 


0.3 


0.61 


0.04 


0.24 


0.03 


NGC 0720 


239 


5 


3.7 


1.8 


0.55 


0.09 


0.40 


0.08 


NGC 0821 


189 


3 


7.3 


1.5 


0.33 


0.02 


0.22 


0.03 


NGC 3608 


178 


3 


6.1 


1.5 


0.38 


0.04 


0.24 


0.04 


NGC 4374 


282 


3 


11.1 


1.5 


0.24 


0.02 


0.28 


0.02 


NGC 4472 


279 


4 


7.8 


1.5 


0.36 


0.04 


0.29 


0.03 


NGC 4552 


252 


3 


10.5 


1.5 


0.36 


0.03 


0.31 


0.02 


NGC 4649 


310 


3 


11.9 


1.5 


0.37 


0.03 


0.32 


0.02 


NGC 4697 


162 


4 


7.1 


1.8 


0.19 


0.04 


0.18 


0.03 


NGC 5638 


154 


3 


7.8 


1.5 


0.32 


0.03 


0.26 


0.03 


NGC 5812 


200 


3 


5.0 


1.1 


0.47 


0.03 


0.26 


0.03 


NGC 5813 


205 


3 


14.9 


2.3 


0.07 


0.04 


0.29 


0.03 


NGC 5831 


161 


3 


2.7 


0.2 


0.61 


0.04 


0.20 


0.03 


NGC 5846 


224 


4 


12.2 


2.4 


0.25 


0.04 


0.29 


0.03 


NGC 6703 


183 


3 


3.9 


1.0 


0.39 


0.05 


0.21 


0.03 



Other galaxies from G93, TMB models: 



NGC 0221 


72 


3 


2.8 


0.7 


0.10 


0.05 


-0.07 


0.01 


NGC 0315 


321 


4 


5.0 


1.5 


0.44 


0.06 


0.32 


0.02 


NGC 0507 


262 


6 


6.9 


2.8 


0.29 


0.07 


0.27 


0.03 


NGC 0547 


236 


4 


8.3 


2.4 


0.34 


0.05 


0.33 


0.02 


NGC 0636 


160 


3 


3.8 


0.7 


0.44 


0.07 


0.18 


0.02 


NGC 1453 


286 


4 


7.1 


1.9 


0.42 


0.06 


0.29 


0.02 


NGC 1600 


315 


4 


7.6 


2.2 


0.47 


0.06 


0.30 


0.02 


NGC 1700 


227 


3 


2.3 


0.7 


0.63 


0.10 


0.17 


0.03 


NGC 2300 


252 


3 


5.5 


1.5 


0.48 


0.05 


0.32 


0.02 


NGC 2778 


154 


3 


5.0 


1.8 


0.40 


0.09 


0.30 


0.03 


NGC 3377 


108 


3 


3.5 


0.8 


0.30 


0.06 


0.27 


0.02 


NGC 3379 


203 


3 


8.0 


1.4 


0.32 


0.03 


0.28 


0.01 


NGC 3818 


173 


4 


5.2 


1.8 


0.47 


0.08 


0.30 


0.03 


NGC 4261 


288 


3 


14.5 


3.3 


0.29 


0.04 


0.27 


0.01 


NGC 4478 


128 


2 


4.3 


2.3 


0.40 


0.10 


0.22 


0.03 


NGC 4489 


47 


4 


2.3 


0.4 


0.24 


0.06 


0.04 


0.02 


NGC 6127 


239 


4 


10.8 


2.2 


0.28 


0.04 


0.30 


0.02 


NGC 6702 


174 


3 


1.4 


0.1 


0.80 


0.07 


0.16 


0.03 


NGC 7052 


274 


4 


11.7 


3.1 


0.27 


0.05 


0.31 


0.02 


NGC 7454 


106 


3 


4.7 


1.0 


0.04 


0.04 


0.13 


0.02 


NGC 7562 


248 


3 


7.1 


1.6 


0.31 


0.04 


0.24 


0.01 


NGC 7619 


300 


3 


13.5 


2.2 


0.31 


0.03 


0.25 


0.01 


NGC 7626 


253 


3 


12.0 


2.4 


0.27 


0.03 


0.32 


0.01 


NGC 7785 


240 


3 


7.9 


2.3 


0.31 


0.05 


0.23 


0.02 
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TABLE 5 
Supplementary Parameters 



Galaxy 


M B 


logr e (pc) 


(v/a), 


7 








Gal 


axics from this work: 






NGC 


0596 


-19.58 


3.41 


0.67 


0.16 




NGC 


1052 


-20.11 


3.53 


1.00 


0.18 




NGC 


1172 


-19.48 


3.68 




-0.01 




NGC 


1209 


1 Q 79 
-iy. i z 


3.40 








NGC 


1400 


-20.16 


3.62 




-0.10 




NGC 


1407 


-21.21 


o.yu 








NGC 


2768 


-20.26 


3.65 




0.24 


-0.054 


NGC 


3115 


-20.30 


3.29 


1.25 


0.52 


-0.069 


NGC 


3193 


-20 28 


3.48 


0.80 


n m 

U.Ul 


-u.uou 


NGC 


3226 


-20.16 


3.83 




0.00 




NGC 


3607 


-21.24 


3.86 


0.92 


0.26 


-0.067 


NGC 


3610 


-20.28 


3.19 


1.10 


0.76 


-0.075 


NGC 


3613 


-19.99 


3.44 


0.84 


0.04 


-0.074 


NGC 


3640 


-19.99 


3.63 


1.48 


-0.10 


-0.043 


NGC 


4168 


-20.18 


3.77 


0.22 


0.17 




NGC 


4365 


-20.34 


3.64 


0.08 


0.07 


-0.058 


NGC 


4473 


-19.77 


3.28 


0.40 


-0.07 


-0.063 


NGC 


4486 


-21.46 


3.90 


0.11 


0.27 


-0.063 


NGC 


4621 


-20.33 


3.55 


0.74 


0.75 


-0.081 


NGC 


5576 


-19.84 


3.29 


0.22 


0.01 


-0.076 


NGC 


5846 


-21.15 


4.02 


0.10 




-0.051 




Galaxies from G93, volume-limited sample: 




NGC 


0584 


-20.39 


3.48 


1.55 


-0.01 




NGC 


0720 


-21.16 


3.75 


0.32 


0.06 




NGC 


0821 


-20.27 


3.64 


0.70 


0.10 




NGC 


3608 


-20.19 


3.61 


0.27 


0.09 


-0.064 


NGC 


4374 


-21.39 


3.68 


0.09 


0.15 


-0.044 


NGC 


4472 


-21.81 


3.93 


0.43 


0.01 


-0.040 


NGC 


4552 


-20.44 


3.36 


0.28 


-0.10 


-0.080 


NGC 


4649 


-21.51 


3.80 


0.42 


0.16 


-0.053 


NGC 


4697 


-20.36 


3.65 


0.71 


0.22 




NGC 


5638 


-20.12 


3.65 


0.73 






NGC 


5812 


-20.40 


3.47 


0.52 


0.59 




NGC 


5813 


-21.20 


3.90 


0.51 


-0.10 


-0.045 


NGC 


5831 


-19.94 


3.57 


0.19 


0.33 


-0.078 


NGC 


6703 


-20.21 


3.50 


0.30 










Other galaxies from G93: 






NGC 


0221 


-15.91 


2.20 


0.89 


0.00 




NGC 


0315 


-22.15 


4.23 


0.09 






NGC 


0507 


-21.74 


4.35 


0.09 


0.00 




NGC 


0547 


-21.20 


3.91 


0.24 






NGC 


0636 


-20.23 


3.45 


1.04 






NGC 


1453 


-21.33 


3.85 


0.62 






NGC 


1600 


-22.23 


4.17 


0.03 


-0.03 




NGC 


2300 


-20.38 


3.65 


0.08 


0.07 




NGC 


2778 


-18.67 


3.34 


0.74 


0.33 




NGC 


3377 


-19.26 


3.28 


0.86 


0.03 


-0.075 


NGC 


3379 


-20.02 


3.27 


0.72 


0.18 


-0.043 


NGC 


3818 


-20.41 


3.58 


0.93 






NGC 


4261 


-21.22 


3.79 


0.10 


0.16 


-0.089 


NGC 


4478 


-19.16 


3.11 


0.84 


-0.10 


-0.103 


NGC 


4489 


-18.46 


3.46 


1.49 






NGC 


6127 


-21.03 


3.82 


0.11 






NGC 


6702 


-20.55 


3.87 


0.18 






NGC 


7052 


-21.14 


3.96 


0.34 


0.16 




NGC 


7454 


-19.34 


3.49 


0.13 






NGC 


7562 


-21.50 


3.86 


0.06 






NGC 


7619 


-21.77 


3.93 


0.53 


-0.02 




NGC 


7626 


-21.03 


3.88 


0.12 


0.47 




NGC 


7785 


-20.91 


3.78 


0.47 


-0.10 





TABLE 6 

Principal Component Analysis Table 



Input PaxcimGters 


% of variance 


Principal Component 


M v , 7, log(t) 


52 


PCI = -0.305Mv - O.68I7 + 0.6661og(t) 


M V , 7, [Z/Hj 


54 


PCI = 0.170M V + 0.7317 + 0.661[Z/H] 


M v , 7, [<*/FeJ 


45 


PCI = -0.258Mv - 0.7027 + 0.663[a:/FeJ 




38 


PC2 = 0.6657 + 0.741 [a/Fe] 


M v , log(r(,), log(t) 


64 


PCI = 0.634Mv - 0.6181og(r 6 ) - 0.4661og(t) 


Mv, log(r ft ), [Z/Hj 


49 


PCI = 0.433Mv - 0.898Iog(r 6 ) 




39 


PC2 = O.lllMv + 0.993[Z/H] 


Mv, log(r ft ), [a/Fe] 


56 


PCI = 0.328Mv - 0.7351og(rs,) - 0.594[a/Fc] 


Mv, 7, log(r 6 ), log(cr), 
log(t), [Z/H], [a/Fe] 


38 


PCI = -0.275Mv - 0.3467 + 0.4601og(r b )+ 
+0.4671og(cr) + 0.4461og(t) - 0.200[Z/H] + 
+0.369[a/Fc] 



25 PC2 = 0.2877 + 0.1601og(r b ) + 0.3731og(cr)- 

-0.3601og(t) + 0.702[Z/H] + 0.351[a/Fe] 



20 PC3 = -0.5917 + 0.5271og(r b ) - 0.3041og(o-)- 

-0.4361og(t) + 0.175[Z/H] - 0.242[a/Fe] 
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Fig. 1. — Spectral indices plotted against galaxy velocity dispersion. The Mg-a relation is clearly apparent, as are similar but 
noiser relations in Ca4227, C24668, and the iron indices. A relation in the opposite direction (smaller index value at larger a) 
is seen in H/3 and H7. 
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Fig. 2. — Comparisons of standard star primary index measurements to the Lick/IDS sample of Wort hev et alJ 1119941) . Different 
point types represent different observing runs. The large solid circle represents HD 51440, an internal standard star for the 
Lick/IDS system and as such the most accurately measured data point in this figure. No calibration correction has been applied 
to these data points. The calibration uncertainty shown here is a major limitation on the overall index uncertainty, as the 
calibration error often dominates the error budget. 
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Fig. 3. — Comparisons of standard star primary index measurements to the Lick/IDS sample of Wort hev et alJ 1119941) . Different 
point types represent different observing runs as above. A constant calibration offset has been applied to each observing run; 
the 2002 November and 2003 January runs have been combined for this purpose due to the limited number of standard stars 
observed in each. No trend with index strength is seen, confirming that a constant offset is sufficient to calibrate the data onto 
the Lick/IDS system. 
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Fig, 4. — Comparisons of uncorrected standard star primary index measurements to the Lick/IDS sample of Worthev et al. 
(1994). The solid point is HD 51440, one of the internal standard stars for the IDS system; as such it is more accurately 
measured than any other star in this data set. Also shown is the typical error in the mean for the index offsets within a given 
observing run. The run to run variations in the mean index values are apparent. 
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Fig. 5. — As in Fig. 2] after correcting the observations to the Lick/IDS system. The offset for each run was determined 
separately except for Runs 4 and 5 which were combined for calibration purposes. 
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Fig. 6. — Stars which were observed repeatedly during the observing program are shown. Unlike in Figs. HI and individual 
index measurements are compared to the average index measurement for each star, not to the Lick/IDS system. Thus this 
shows the internal error in each index, independent of the uncertainty in calibrating to the Lick/IDS system. 
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Fig. 7. — The calibrated primary index measurements for all standard stars are plotted against the standard value for that 
star on the Lick/IDS system. Open circles are ordinary standard stars from this work, and the solid circle as above is the 
Lick/IDS internal standard star HD 51440. Several stars were observed as part of both G93 and this study; the calibrated G93 
measurements for those stars are shown as open triangles. The adopted calibration is in good agreement with both the Lick/IDS 
system and the G93 calibration. 
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Fig. 8. — Comparisons of measured velocity dispersions between galaxies in this s tudy and in other d ata sets. Open tria ngles 
are compared with G 93. Black open circles, x's, and solid circles are compared with lFaber et all (ll989T ). lFaber et alJ ljl997f) . and 
iDenicolo et, al.l (f2005Tl respectively. The velocity dispersion measurements from this study are in good agreement with previous 
work. 
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Fig. 9. — Comparisons of measured spectral indices between gal axies in this study an d in other data sets. Points in large type 
are compared with G93. Points in small type are compared with iDenicolo et alJ (|2fl05). 
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Fig. 10. — Comparisons of measured vel o city d ispersion and spectral indices between galaxies in G93 and in IDenicolo et all 
J200fJl. The fact tha t the IDenicolo et"all j200fifl data set deviates from the G93 in each quantity in the same way that the 
IDenicolo et, al.l IMPORT) data set deviates from the volume-limited data set (Fig. El lends confidence that the volume-limited data 
set has been successfully calibrated onto the Lick/IDS system. 
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[MgFe]' [A] 

Fig. 11. — The galaxies in the volume-limited sample are plotted on age-metallicity grids from TMB. The horizontal axis is the 
[MgFe]' index defined by TMB to be independent of [a/Fe]. Galaxies in normal type are from this study, while galaxies in italics 
are members of the volume- limited sample from the G93 data set. A dashed line connects the measurements of NGC 5846 from 
the two data sets; this galaxy provides a direct empirical test of the accuracy of the index measurements and consistency of the 
Lick/IDS calibration. Although the H/3 measurements are in excellent agreement between the two studies, the Mgfo, Fe5270, 
and Fe5335 index measurements are discrepant by more than two standard deviations. Much of this discrepancy results from 
the unusually large difference in measured velocity dispersion between the two studies. Typical error bars are shown in the 
lower left. The volume- limited sample is in agreement with previous work, showing a wide range of ages and a smaller range of 
supersolar metallicities. 
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Fig. 12. — The galaxies in the volume-limited sample are plotted against iron and magnesium lines to measure [a/Fe]. Models 
from TMB are shown, at 3 Gyr age (solid lines) and 8 Gyr age (dotted lines), and [a/Fe] values of 0.0 and +0.3 dex. Galaxies 
are presented as in 

Fig. El Typical error bars are shown in the lower right corner. The two sets of models shown illustrate the 
relatively small effect of age on [a/Fe] measurements. 
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Fig. 13. — Derived physical parameters for the G 93 sample are compared. Parameters from the TMB models are on the vertical 
axes, while parameters from lTrager et all (|2Q0Qa) are on the horizontal axes. Age is compared in the upper left, metallicity in 
the lower left, and [a/Fe] in the lower right. Representative error bars are shown. The three galaxies with similar [a/Fe] using 
both models are NGC 1700, NGC 5831, and NGC 584, the three youngest galaxies plotted. The qualitative difference between 
the two models is negligible. TMB models result in a constant offset toward larger [a/Fe], which in turn results in larger [Z/H]. 
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Fig. 14. — Comparison of H/3 and H7F index measurements. The galaxies in small type have large emission corrections and 
correspondingly large uncertainties. Typical uncertainties for the rest of the galaxy sample are presented in the upper left. With 
the exception of four outliers, the Balmer index measurements follow a tight linear relation. 
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Fig. 15. — The galaxy sample from this work is plotted on the H/3 vs. (Fe) plane, al ong with two set s of model grids. Models 
in solid type are from TMB with [a/Fe]= 0.0, while models in dashed type are from ISchiavon I (2005). In both sets of models 
only the most metal-rich isometallicity line (TMB: [Z/H]= +0.67; Schiavon: [Fc/H]= +0.2) is plotted (dotted lines). Isochrone 
lines for each set of models are labeled. As above, galaxies in small type have large emission corrections. Typical uncertainties 
for the remaining galaxies are presented in the lower left. A consistent age offset is apparent between the two sets of models, in 
the sense that the Schiavon models measure older ages by ~ 2 Gyr compared to the TMB models. 
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Fig. 16. — Galaxies and dashed models are as in Fig. 1151 Models in solid type are from IThomas. Maraston. fc Kornl |2004) 
with [a/Fe]= 0.0. The two sets of models are substantially different from each other in the H7F vs. (Fe) plane. This is primarily 
due to the IThomas et al] 1120041) models predicting larger H7F equivalent widths by ~ 0.4 A compared to the Schiavon models . 
As a result, age measurements using the Schiavon models are younger by several Gyr compared to the IThomas et al] i2004f) 
models. 
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Fig. 17. — Ages derived using the model grids in Figs. fTol and 1161 are compared. Measurements in the left panel are based 
on the Thomas et al. models, while measurements in the right panel are based on the Schiavon models. As above, galaxies 
in small type have large emission corrections. Typical age uncertainties for the remaining galaxies are presented in the lower 
right. Dotted lines are drawn at Age = 14 Gyr on each axis. The galaxies with inconsistent Balmer measurements (outliers 
in Fig. 1140 result in discrepant ages using both sets of models. The two sets of models are in qualitative agreement, though 
as expected from Figs. [T51 and 1161 thev are in systematic quantitative disagreement, typically by 2-4 Gyr. The Thomas et al. 
models result in older ages on the H/3 grids (Fig. 1150 and younger ages on the H7F grids (Fig. 1160 . Neither set of models gives 
fully consistent ages between the two grids. The Schiavon models are more consistant than the Thomas et al. models at young 
and intermediate ages, while the Thomas et al. models are more consistent for extremely old (Age > 12 Gyr) galaxies. For the 
majority of galaxies, the difference between ages derived using different sets of models is significantly greater than the difference 
between ages derived using different Balmer indices. As a result, the conclusion is that the choice of model is more important 
than the choice of age-sensitive Balmer index in measuring SSP parameters. 
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Fig. 18. — Principal component projection of the metallicity hyper-plane of TFWG. Solid circles are from this study. The 
remaining points are from G93; open triangles are in the volume-limited sample while skeletal symbols are not. Error ellipses 
are shown for the measured index and calibration errors for this study (solid ellipses). The projection of each SSP parameter 
in principal component space is shown. PCI is primarily sensitive to log a and [a/Fe], with some metallicity sensitivity. PC2 
is sensitive to age and metallicity, while PC3 is depends weakly on both log a and [a/Fe]. Conceptually, one can think of PCI 
measuring a galaxy's position along the Mg-cr relation, PC2 measuring position along an age-met allicity relation, and PC3 
measuring deviations from the Mg-a relation. The volume-limited and G93 samples are consistent with being drawn from the 
same distributions in each principal component axis. 
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[Z/H] 

Fig. 19. — The Z-plane as defined in TFWG (solid line). Points are as in Fig 1181 A typical la error ellipse is shown in the 
upper left. The normalization of this relation between velocity dispersion, age, and metallicity depends on the models used, as 
would be expected based on Fig. 1131 The difference in models readily expl ains the offset betwe en the TFWG fit and the data. 
The most metal-rich galaxy is NGC 3610, which is examined in detail in Howcll ~t alJ i2004f) . The volume-limited and G93 
samples are consistent with being drawn from the same distribution along this edge-on projection of the Z-plane. 
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Fig. 20. — Comparison between parameters relating to the average stellar composition (age, [Z/H], [a/Fe]) and parameters 
relating to galaxy size (a, Mb, log r e ), shape ((n/cr)*), density profile (the power-law slope 7), and color gradient. Points are as 
in Fig ll8l Several well-known trends are apparent. The [a/Fe]-a relation is directly analogous to the Mg-cr relation (Fig. |Q. It 
is also clear that the faintest and least massive galaxies are all relatively young, while brighter and more massive galaxies span 
a wide range of ages. 
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Fig. 21. — The face-on projections of the Fundamental Plane (top panel) and Z-plane (bottom panel) are shown. Galaxies in 
normal type are from this work, while those from G93 are in italic type. No mapping from one plane to the other is found. 



